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Abstract The atmosphere of the Sun is characterized by a complex interplay of compet- 
ing physical processes: convection, radiation, conduction, and magnetic fields. The most 
obvious imprint of the solar convection and its overshooting in the low atmosphere is the 
granulation pattern. Beside this dominating scale there is a more or less smooth distribu- 
tion of spatial scales, both towards smaller and larger scales, making the Sun essentially a 
multi-scale object. Convection and overshooting give the photosphere its face but also act 
as drivers for the layers above, namely the chromosphere and corona. The magnetic field 
configuration effectively couples the atmospheric layers on a multitude of spatial scales, 
for instance in the form of loops that are anchored in the convection zone and continue 
through the atmosphere up into the chromosphere and corona. The magnetic field is also an 
important structuring agent for the small, granulation-size scales, although (hydrodynamic) 
shock waves also play an important role — especially in the internetwork atmosphere where 
mostly weak fields prevail. Based on recent results from observations and numerical simu- 
lations, we attempt to present a comprehensive picture of the atmosphere of the quiet Sun as 
a highly intermittent and dynamic system. 
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1 Introduction 

Observations of the solar atmosphere reveal a wealth of different phenomena, which oc- 
cur over an extended range of different temporal and spatial scales. This is not surprising, 
considering the fact that already basic parameters such as gas density and temperature span 
many orders of magnitude, from the convection zone below the photosphere to the corona. 
At a first look, it may thus appear rather hopeless to construct an overall picture that can 
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account for all the phenomena. At a closer look, however, many connections between ap- 
parently independent phenomena can be found, ultimately implying a multitude of cou- 
plings through the atmosphere. In addition, there seems to be a hierarchical arrangement 
of approximately selfsimilar convective motions, with the granulation pattern embedded in 
increasingly larger meso- and supergranulation patterns. 

The key to a comprehensive picture of the solar atmosphere thus lies in relaxing too strict 
and oversimplified concepts, even when they are didactically nicer than the reality. The solar 
atmosphere should not be seen as a static stack of layers but rather as intermittent domains 
that are dynamically coupled together. One example is magnetic flux structures (or "flux 
tubes") fanning out with a wine-glass geometry. Such regular building blocks put certain 
constraints on the implied atmospheric structure, which can make it difficult to fit in other 
observational findings. Accepting that magnetic field structures are far less regular offers 
room for a more generally valid comprehensive picture. This trend became more and more 
obvious during the recent years, both from the observational and theoretical side (see, e.g., 
|Carlsson|2007l|Gudiksen|200^|Hansteen|2007[|E^ and 
many more). 

The advantages of a relaxed picture can be seen from the example of the quiet Sun 
chromosphere above internetwork regions, which in itself is a complex and intriguing phe- 
nomenon (see, e.g., |Judge|200^|Rutten|2^0^[DePbntieu et al.|2004]|Lites et al.|1999| and 
many more). Despite tremendous progress, there are still many open questions concerning 
its structure, dynamics and energy balance. Recent observations now prove - beyond any 
doubt - the chromosphere to be a highly dynamic and intermittent layer. The internetwork 
chromosphere is the product of a dynamic interplay of shock waves and magnetic fields. This 
picture, which was already suggested by many earlier investigations, offers a key to resolve 
some apparent contradictions that lead to much confusion in the past. A prominent example 
concerns the observation of carbon monoxide (see, e.g., Ayres 2002J, which now can be ex- 
plained as an integral part of a dynamic and intermittent atmosphere (fWedemeyer-Bo hm &| 
Steffen 2007; Wedemeyer-Bohm et al. 2006, 2005a). And still the chromosphere cannot be 
investigated without also taking into account the layers above and below. The shock waves, 
which are so essential at least for the lowest, weak-field parts of the chromosphere, are gen- 
erated in the layers below, while significant amounts of mass and energy are exchanged 
between the chromosphere and the corona above. Obviously, the whole atmosphere must be 
seen as an integral phenomenon. 

In the following sections, we report on a selection of results from observations and nu- 
merical simulations, which will help us put together an updated, revised view of the structure 
of the quiet Sun atmosphere. 



2 The Sun - a multi-scale object 

An overarching point in this discussion is the fact that the Sun is fundamentally a multi- 
scale object. This is a major difficulty for modeling and understanding, since it requires 
(computationally expensive) modeling over a large range of scales. 

But the Sun also displays aspects of self-similarity and scale invariance in several re- 
spects, which on the other hand helps a lot. To illustrate the self-similarity, Fig. [T^ shows 
temperature patterns in horizontal planes in a large scale simulation of solar convection 
(Zhao et al. 2007), and Fig. [TJ> shows patterns of vertical velocities from the same sim- 
ulation. The temperature patterns show very intermittent cold structures, embedded in a 
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Fig. 1 Large scale solar convection (48 Mm X 48 Mm X 20 Mm), a) Temperature (top four panels) and b) ve- 
locity (bottom four panels) patterns at four different depths. Temperature is shown on a linear scale. Velocities 
are rendered with positive (downward) values blue and negative (upward) velocities yellow. A narrow band 
near zero velocity is rendered in grey. 



Fig. 2 A side view of a 48 Mm X 48 Mm X 20 Mm simulation, showing velocity streamlines, with brightness 
increasing with increasing magnitude. Up- and down-flows are rendered in blue and red, respectively. 

background of horizontally nearly constant temperature (images of entropy would look es- 
sentially identical, with near-constancy also in the vertical direction). The set of panels also 
shows that the pattern scales increase systematically with depth. 

Figure [lj>, on the other hand, which displays vertical velocity on a color scale that 
changes from yellow to blue with sign (with a narrow band of grey for velocities near zero) 
gives a completely different impression. With this rendering choice one can see that, at least 
from a morphological point of view, the patterns at different depths are quite similar. Dis- 
playing in this way, signed velocity reveals that the sharply defined dark (cold) patterns 
in Fig. [T^ indeed correspond to the strongest downward velocities, but that there are also 
relatively broad areas of much milder downflows. This shows that, as the ascending gas is 
forced (by mass conservation) to overturn, it does so at first gently, then to finally be ac- 
celerated more strongly by the positive feedback that comes from merging with colder gas 
from above. 

At the visible surface the horizontal velocity patterns from various depths are super- 
posed. This happens because the depth dependencies of the large scale horizontal velocity 
patterns are rather weak; at least over depth intervals small compared to their horizontal 
extents. On the other hand, as illustrated by Fig.[T] the dominant scale becomes smaller for 
layers increasingly close to the surface. As the amplitudes of these smaller scales are larger, 
they mask the presence of the larger scale patterns, whose presence, however, can still be 
revealed, e.g., with Fourier analysis or with low-pass filtering. The hierarchy is illustrated in 
a side view in Fig. [2] The combination of streamlines and colors illustrate how near-surface, 
small scales fluctuations are carried along in larger scale flows. 

The hierarchy of scales displayed reveals no particular preferred scale above the gran- 
ular one; the transition to larger and larger scales with depth is smooth. A direct way to 
illustrate this from observations is to use power spectra of solar velocities, as observed with 
SOHO/MID (Georgobiani et al. 2007). Figure [3] shows velocity (mainly horizontal) as a 
function of size, produced by filtering the velocity power observed by MDI into sonic and 
sub-sonic parts. The velocity spectrum displayed is produced by then taking the square root 
of the velocity power times wave number; this is a quantity — a velocity spectrum — that 
nicely illustrates the dependence of velocity amplitude on size. Note that there is very little 
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Fig. 3 Solar velocity spectrum: A subsonic (7 km s~ ) filter has been used to separate the velocity into oscil- 
latory (grey) and convective (black) components. The dashed line shows the convective component resulting 
from first taking a 24h average. Adapted from Georgobiani et al. 1 2007 1. 




Fig. 4 Solar horizontal velocities observed with SOHO/MDI. A patch some distance away from solar center 
has been compensated for projection effects and filtered to effective resolutions that differ by factors of 2. 
Which is which? 

(less than a factor two) extra power at scales traditionally associated with supergranula- 
tion, and that there is a smooth and increasing distribution of velocity amplitude across the 
"meso-granulation" range to granulation scales. The same behavior is found in large scale 
numerical simulations ( Georgo biani et al.[2 007 ). 

Supergranulation patterns can be brought forward by averaging over either time or 
space; the dashed line in Fig. [3] shows the effect of a 24-hour time average. A low-pass 
wavenumber spatial filter has a similar effect; it cuts away the larger amplitudes at smaller 
scales and exposes aspects of the underlying larger scale pattern. The relatively distinct ap- 
pearance of a supergranulation scale network in magnetically related diagnostics indicates 
that the transport and diffusion of magnetic field structures at the solar surface results in 
what is effectively a low-pass wavenumber filter. 
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As shown by Fig. [4] it is practically impossible to tell the difference between velocity 
patterns on different scales, once they are filtered to effectively the same resolution. As 
illustrated by Stenfio & Holzreuter ( 2002 2003 1 magnetic field patterns and distributions 
also show a degree of self- similarity. 

The magnetic energy equation 

= -V-(ExB)-u-(jxB)-e Joule (1) 

illustrates that balance of the magnetic energy at each depth is achieved by Lorentz force 
work (by the flow on the field) being used to balance magnetic dissipation, with net magnetic 
energy transported up or down by the Poynting flux, E x B. As shown by the work of 
Vogler & Schiissler (2007) and St einer et al.| l |2008| > the actual direction of net transport is 
systematically downwards, at least below the solar surface. It appears likely that there is net 
dynamo action at each depth in the convection zone, with net magnetic energy delivered 
to the next layer down. This naturally leads to a pile up near the bottom of the convection 
zone. The downward transport, which is known from direct studies ( jTobias et al.|1998[|2001| 
|Dorch &"N ordlund 2001 ), is often referred to as "turbulent pumping", and is associated with 
the asymmetry between up and downflows (illustrated in Fig.fTJ)). 

On the largest scales (largest depths), and only there, differential rotation enables a large 
scale global dynamo action, with patterns clearly controlled by being stretched out by dif- 
ferential rotation. Buoyancy eventually pushes the fields back up. 

Another evidence for self-similarity comes from the power law behavior of flare energy 
distribution as a function of time. This behavior is also recovered in numerical simulations 
of 3D magnetic reconnection (Galsgaard & Nordlund 1996). These are signs that magnetic 
reconnection occurs in a multi-scale hierarchy, where magnetic dissipation at large magnetic 
Reynolds number (low resistivity) creates a hierarchy. Large scale structures generate sub- 
sidiary small scale structures, which do it again (on shorter time scales) and again, until the 
spatial scales are small enough to support the dissipation. 

Note the remarkable and wonderful argument, made already by Parker a long time ago, 
which shows that driven magnetic dissipation must, if anything, increase with decreasing 
resistivity — quite contrary to naive expectations. This has been verified in numerical exper- 
iments by at least three different groups (Galsgaard & Nordlund 1996; Hen drix et al.| 1996, 
IDmitruk & G6mez|1999| l. 

The chromosphere and corona are likely to be heated in much the same way, as is illus- 
trated by the well known flux-flux relations between coronal and chromospheric diagnostic. 
It is hard to even avoid, as in models of coronal heating there is a tendency of dumping much 
more energy in the chromosphere as a side effect (Gudiksen & Nordlund 2002, 2 005b) >. 

As pointed out by Phil Judge: The chromosphere is not a mess; the upper chromosphere 
looks nearly force-free like the corona, whereas the lower chromosphere is less force-free. 
Complexity comes from both the temperature and density. A central question is: What drives 
the flows (particularly the cool upflows)? 

Semi-realistic models of coupling of the horizontal photospheric velocity field to the 
corona were first computed by Gudiksen & Nordlund ( 2005b ), who showed that a correctly 
normalized photospheric (model) velocity field injects sufficient power into the corona to 
create and maintain coronal temperatures (cf. Fig. [5}. 

The mechanism is, in the absence of explicit flux emergence, essentially the 'braiding 
mechanism' introduced by Parker ( 1972 1981 , 1983). The heating is quite intermittent and 
drives up- and down-flows along magnetic field lines into and out of the corona. Emulated 




7 



Temperature 




igffi|fc=3 B: 4.5 Mm 
P^=^ A: 0.3 Mm 



Fig. 5 The panels show the vertical magnetic field (far left), the ratio (often referred to with the symbol a) 
of the magnitude of the electric current along the magnetic field to the magnitude of the magnetic field itself 
(left), the gas pressure (right) and the log. temperature (far right). The positions of the cutting planes are 
indicated in the inset at bottom right. Adapted from Gudiksen & Nordlund ( 2005b I. 




Fig. 6 The Si IV 1 394 A Doppler shift (left) and emission measure (mid) as they would be observed from 
above, and the emission measure projected along the Y-axis (right). Positive (downward) Doppler shifts are 
in red and negative Doppler shifts are in blue. Adapted from Peter et al. (2 004) . 



TRACE images and animations show qualitative agreement with observations l |Peter et aT] 
2004). Silicon IV, for example, picks up cooling condensations (cf. Fig. |6f . In addition, 
spectral lines formed at different temperatures show semi-quantitative agreement of the de- 
pendency with depth of their Doppler shifts and mean emission measures (cf. Figs. 7 and 
9, |Peter et al.|20 06). The differential emission measure is a 'fingerprint' type diagnostic, in 
much the same way as spectral line asymmetries are for photospheric spectral lines ( |Dravins| 
|et al.|1981||Dravins & Nordlund|1990||Asplund et al.|2000a|b|>. Subsequently ther e has been 
much progress due to the work of the Oslo group ( |Hansteen & Gudiksen|2005[|Hansteen| 
|et al.|2006|[2007{|Martmez-Sykora et al.|2008| > - cf. also the discussion in Sect.|4.3| 
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3 Observations - Measuring the magnetic field in the solar atmosphere 

The Ho- line core images in Figs.[7]i and [8] show a well-known but still barely understood 
and intricate picture: fibrils that spread from regions of enhanced magnetic field strength, 
occasionally connecting to neighboring regions or apparently fading in between (see e.g. 
|Rutten|2007| l. The structure gradually changes as one goes from line center into the wings, 
as the corresponding intensity is due to lower layers. Finally, a mesh-like background pat- 
tern shines through in the internetwork regions. It is most likely due to reversed granulation 
in the middle photosphere with some possible contributions from the low chromosphere. 
The gradual change of the pattern in Ho- with wavelength gives some clues about the atmo- 
spheric structure, in particular the magnetic field in the chromosphere (the"canopy" field), 
and definitely shows us that the photosphere and chromosphere are coupled via magnetic 
fields on medium to large spatial scales and via fields and shock waves on the small scales. 

Therefore, the understanding of the coupling between photosphere and corona is inti- 
mately connected to the measurement of the chromospheric magnetic field. The following 
subsections exemplify the difficulties of chromospheric magnetic field measurements and 
present promising approaches to determine the vector magnetic field of the chromosphere. 



3.1 Improving magnetic field extrapolations 

Measuring the magnetic field in the photosphere has a long tradition [Hale| (|T908 1. After 100 
years of solar magnetic field measurements, the level of sophistication, both in terms of in- 
strumentation and in analysis technique, has reached a very high level of maturity. Numerous 
telescopes, ground based and space born (e.g. GONG, MDI), investigate the global struc- 
ture of the solar magnetic field on a routine basis. High resolution measurements allow the 
characterization of magnetic elements with a size in the 100 km range (e.g., SST - CRISP, 
Hinode SP). With HMI on SDO and the balloon-borne 1 m telescope Sunrise (launch: sum- 
mer 2009) a major improvement in the determination of photospheric magnetic fields will be 
achieved in both directions - the global magnetic field configuration as well as the smallest 
scale structures down to a size of 25 km. 

The availability of such high quality photospheric vector magnetograms and the low 
plasma-/? in the chromosphere are the basic ingredients needed for reliable, force-free mag- 
netic field extrapolations. Starting with Sakurai ( 1981 ) these extrapolations nowadays have 
reached a high level of sophistication (see reviews by Sakurai 1989, [Amari et aL||1997| 
Wiegelmann 2008). To further improve the accuracy of the chromospheric magnetic field 
extrapolations additional information on the complex structure of the chromosphere must 
be taken into account. One of the most promising advances in this direction was proposed 
by |Wiegelmann et al.| < [2008] >: the basic assumption for applying non-linear, force-free mag- 
netic field extrapolations is the force-freeness of the photospheric vector magnetograms. 
Measured magnetograms do not fulfill this requirement, therefore a preprocessing of the 
measured data is required. Wiegelmann et al. (2008 ) developed a minimization procedure 
that yields a more chromosphere-like field by including the field direction information con- 
tained in, e.g., chromospheric Ha images. Including this information into the extrapolation 
algorithm significantly enhances the reliability of the extrapolations. 
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Fig. 7 Observations of a quiet Sun region close to disc-centre: a) Ca II H wide band, b) Ha wide band image 
(FWHM 0.8 nm), c) Ca II H wing (396.5 nm), d) Ca II H inner wing (close to line core), e) Ca II H core, 
f) Ha blue wing at -35 pm, g) Ha red wing +35 pm, h) Ha line core. The observations were carried out with 
the Swedish 1-m Solar Telescope (SST). Data courtesy: L. Rouppe van der Voort (University of Oslo). 
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Fig. 8 Ha fine structure at -800, -600, -400, -200 mA from line center recorded with the Dutch Open Tele- 
scope (DOT) on October 4, 2005 (taken from Rutten 2007 1. Cell -spanning fibrils are visible around line center 
(right). The decreasing line opacity in the blue wing of the line opens the view to the solar photosphere, in- 
tercepted by dark fibrils resulting from the Doppler shift of the line core. 



3.2 Direct measurements of the chromospheric magnetic fields 



Measurement techniques for chromospheric magnetic fields have to overcome a variety of 
hurdles: (i) the plasma density is several orders of magnitude lower than in the photosphere, 
(ii) the energy transport is dominated by radiation, (iii) the magnetic field strength is on 
average lower than in the photosphere, and (iv) anisotropic illumination induces population 
imbalances between atomic sublevels that are modified by weak magnetic fields. The low 
plasma density leads to weak signals in the absorption (on-disk observations) or emission 
(off-limb observations) of spectral lines. The absorption signatures of chromospheric lines 
often show a strong photospheric contribution. Only highly spectrally resolved observations 
of the line core carry the chromospheric information. As a consequence of the low density, 
the simplifying assumption of local thermodynamic equilibrium breaks down. The interpre- 
tation of the observations is thus by far more involved than in the case of photospheric obser- 
vations. Additionally, the low chromospheric magnetic field strengths weakens the Zeeman 
signals in spectral lines. Scattering polarization and its modification by the Hanle effect in- 
troduce an additional complication in the analysis of the polarization signal of spectral lines. 

During the last decade major progress has been achieved in circumventing these hurdles. 
Radio observations are able to determine the magnetic field strength in and around active re- 
gions (see review by Lee 2007 1. Acoustic mapping techniques (Finsterle et al. 2004) use the 
reflection of high-frequency acoustic waves (mHz-range) from the region in the atmosphere 
where the gas pressure and the magnetic pressure are equal to reveal the structure of the 
magnetic canopy. The biggest leap in the direct determination of chromospheric magnetic 
fields was achieved by combining state of the art instrumentation for full Stokes polarime- 
try with recent progress in atomic physics. Bommier ( 1980), Landi DeglTnnocenti ( 1982). 
|Stenflo & Kelie^fl99"7] l and |Trujillo Bueno et al.| ( |2002| i opened a new diagnostic window in 
solar physics: magnetic fields influence the strength and the direction of the linear polariza- 
tion resulting from atomic or scattering polarization. This effect, discovered by Hanle ( 1924 
|(reprint 1991) , allows the determination of the magnetic vector from Milligauss to several 
tens of Gauss, a range not accessible by Zeeman diagnostics. 

The following sections describe examples of measurements in this new diagnostic win- 
dow, focused around two of the most popular spectral lines for combined Hanle and Zeeman 
measurements: the triplet of He i 10830 A and the He i D3 5876 A multiplet. The formation 
of these lines requires ionization of para-He by ultraviolet radiation or collisions, followed 
by recombination to populate the lower sates of ortho-He. Since the main source for the ul- 
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Fig. 9 Magnetic field measurement in the He i D3 line: the magnetic field vector determined by a combined 
Hanle and Zeeman diagnostic traces the visible structures in the Ha- slit-jaw image (bottom, adapted from 
|L6pez Ariste & Casini|2005) . 



traviolet radiation is the corona, these He 1 lines lacks almost any photospheric contribution. 
Additionally, they are (generally) optically thin and narrow, allowing the use of rather sim- 
ple analysis techniques, like Milne-Eddington inversions of the radiative transfer equations 
( [Solanki et al.|2003[|Lagg et al.|2004||Lagg et al.|2007). Wi th the inversion code HAZEL 
(HAnle and ZEeman Light, see |Asensio Ramos et al.||2008] l, involving the joint action of 
atomic level polarization and the Hanle and Zeeman effect in these lines, a standard tool for 
the analysis of Stokes spectra is now available. 



3.2.1 Spicules 

Spicules are an ubiquitous phenomenon on the Sun. At any time, the number of these needle- 
like structures on the Sun is on the order of 4 x 10 5 . These dynamic and short lived features 
(lifetimes typically 5-10 minutes) can be considered as magnetic tunnels through which the 
refueling of the coronal plasma takes place 1 Athay|2000"|l. Hig h cadence Hinode SOT obser- 



vations in Ca 11 H ( Okan ioto et al.|2 007 ; Suematsu et al. 2007 ) revealed details in terms of 
size and dynamics and led to the discovery of a new type of spicules (type II spicules, De 
|Pontieu et al.|2007b) with shorter lifetimes (10-150 s), smaller diameters (<200km com- 
pared to <500km for type I spicules), and faster rise times. According to |De Po ntieu 
et al. (2007c), they (i) act as tracers for Alfven waves with amplitudes of the order of 10 
to 25 kms -1 and (ii) carry, in principle, enough energy to play an important role for heating 
the quiet Sun corona and to accelerate the solar wind. See also Sect.|5.2| 
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Fig. 10 Intensity image in the center of the He i 10830 A line and derived azimuthal and inclination angle of 
the magnetic field (adapted from Merenda et al. 2007 1. 



Measurements of the magnetic field of spicules, both type I and type II, are essential 
for the understanding of this phenomenon. Truji llo Bueno et"aL] l |2005| l were the first to 
directly demonstrate the existence of magnetized, spicular material. Full Stokes po larimetric 
data in the He i 10830 A line, obtained with the Tenerife Infrared Polarimeter ( Collados 



et al. 1999), were analyzed by solving the radiative transfer equation assuming an optically 
thick atmosphere. The application of a combined Hanle and Zeeman diagnostic revealed 
a magnetic field strength for the observed type I spicule of 10 G and an inclination angle 
of 37° at a height of 2000 km above the photosphere. The authors state that 10 G is the 
typical field strengths for spicules at this height, but significantly stronger fields may also 
be present. This result agrees with the measurements from Lopez Ariste & Casini ( 2005 1 
using full Stokes polarimetry in the He i D3 line. They find field strengths not higher than 
40 G and a good correlation between the magnetic field orientation and the visible structure 
in Ha (see Fig. [9]l. An ind ependent confirmation of these measurements was presented by 
|Socas-Navarro & Elmore| ( [2005[ l by using full Stokes observations from SPINOR (Spectro- 
Polarimeter for INfrared and Optical Regions, at the Dunn Solar Telescope). Their multi- 
line approach removes the dependence of the strength of the Hanle signals on the zero-field 
polarization produced by the scattering of anisotropic radiation in the higher atmosphere. 



3.2.2 Prominences and filaments 



The spectacular eruptions of prominences and filaments and the resulting coronal mass ejec- 
tions (CME) can cause sudden changes in the terrestrial magnetosphere. A typical C ME re- 
leases an energy of 10 25 J and 10 12 kg of solar material into the interplanetary space I jHarri- 
|son|1 994). Before eruption, solar magnetic field holds this dense and relatively cool material 
in the hot coronal environment and supports it against the solar gravity for time periods as 
long as weeks. The knowledge of the magnetic field within these structures therefore is of 
great interest to understand the mechanisms leading to a possible eruption. 

|Casini et al.| < |2003j ) were the first to present magnetic maps of prominences using full 
Stokes polarimetry in the He 1 D3 line. Their results confirm previous measurements of the 
average field in prominences, ranging between 10 and 20 G and oriented horizontally with 
respect to the solar surface. However, they also find the presence of organized structures in 
the prominence plasma embedded in magnetic field that are significantly larger than average 
(50 G and higher). Merenda et al. ( 2007 1 extended this work to include the forward scattering 
case and applied it to a filament located at disk center and obtained the first magnetic maps of 
a filament. In this preliminary work they restricted their analysis to the saturated Hanle effect 
regime between 10 and 100 G. Here the linear polarization is only sensitive to the direction 
of the magnetic field and does not change with intensity variations. The results are reliable 
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Fig. 11 Speckle-reconstructed, narrow band image (contrast-enhanced) of a plage region observed in the line 
core of Ca K using the SST jPietarila et al.|2008| cf. |Rouppe van der Voort et al.|2005) . The Ca K fibrils 
extend over quiet Sun regions. The mesh-like background pattern is nevertheless dominated by the reversed 
granulation pattern in the middle photosphere (cf. Fig.[JJ. 



maps for the azimuth and inclination angle for the magnetic field (see Fig . 1 1 0[> . In agreement 
with Casini et al. (2003) they find horizontal fields in the central part of the filament and 
a change of the azimuth according to the orientation of the main axis of the filament. In 
order to detect, for example, the small-scale and rapidly moving filaments mentioned in 
Sects. |4~2| and |5.3| significant improvements in signal to noise ratio and temporal resolution 
of polarimetric observations in this spectral line are required. The complex magnetic field 
and velocity structure of an erupting filament in the He i 10830 A line was analyzed by 
Sass o et aT] < |2007[ >: besides the magnetic field topology they identify the presence of up 
to 5 different atmospheric components, distinguished by their velocities ranging from -50 
to +100 kms -1 , within the resolution element of approximately l."5. This measurement 
clearly demonstrates the fibrilar structure of the chromosphere (see also Lagg et al. 2007 1 
and the need for higher spatial resolution measurements in this line. 

3.2.3 Canopy 

Following previous work by W. Livingston, |Gabriel| ( |1976| l introduced the term canopy to 
explain the emission measures of chromospheric and transition region UV lines. In the "clas- 
sical" picture, the magnetic pressure wins over the gas pressure with increasing height, so 
that the magnetic flux concentrations rooted in the network expand and cover the inter- 
network cells with horizontal fields (see Sect. [5] for an updated view). Giovanell i & Jones| 
( 1982|l; |Jones & Giovanelli| ( [T9 82l performed detailed studies of the magnetic canopy close 
to the limb by determining magnetograms using chromospheric spectral lines like the Ca n 
triplet at 8542 A or the Mg ib2 line at 5173 A. These magnetograms are characterized by a 
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Fig. 12 Measurement of the magnetic field over a supergranular cell in the photosphere and the chromosphere 
(German Vacuum Tower Telescope, Tenerife Infrared Polarimeter 2, May 10, 2008): continuum close to the 
10830 A line (top left), Stokes V signal integrated over the red wing of the photospheric Si i 10827 A line 
(bottom left), Stokes U and V signal integrated over red wing of the chromospheric He 1 1 0830 A line (top and 
bottom right, respectively). The chromospheric maps suggest the presence of magnetic structures organized 
on mesogranular scales within the supergranular cell outlined by the photospheric Stokes V map. 



polarity inversion line parallel to the limb, on either side surrounded by diffuse fields above 
the internetwork region (see |Steiner & M urdin 2000J for a sketch of the magnetic configu- 
ration). 

Especially during the last decade diagnostic tools involving the Hanle effect signifi- 
cantly improved the possibilities to characterize the canopy fields. Using spectropolarimet- 
ric data in the Sr n 4078 A line "Hanle histograms", showing the statistical distributions of 
the Hanle rotation and depolarization effects, Bia nda et al.| lfl998 1 determined the magnetic 
field strength of horizontal, canopy-like fields to be in the range of 5 to 10 G. The first spatial 
mapping of Hanle and Zeeman ( jStenflo et al.|2 002 ) effect revealed details of canopy fields in 
a semi-quiet region measured close to the limb in the Na i D1-D2 system. The authors found 
direct evidence for horizontal magnetic fields, slightly stronger than the field strengths de- 
termined by |Bianda etaT] (19981 (25-35 G), that remain coherent over a spatial scale of at 
least three supergranules. 

The concept of a magnetic canopy around sunspots and in active regions is well estab- 
lished. Over quiet regions, the formation of this layer of horizontal fields is matter of debate: 
|Schrijver & Titl e ( 2003 ) showed that concentrations of magnetic flux in the network in the 
order of a few tens of Mx cm -2 will destroy the classical, wineglass-shaped magnetic field 
topology. Such flux concentrations, suggested by simulations, were identified by Trujillo 
Bueno et al. ( 2004 1 in terms of ubiquitous tangled magnetic field with an average strength 
of ~ 1 30 G, much stronger in the intergranular regions of solar surface convection than in the 
granular regions. A significant fraction of this hidden magnetic flux has now been clearly 
identified with the spectropolarimeter of the Hinode spacecraft (Lites et al. 2008a). How- 
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ever, narrow-band (0. 1 A) observations in the Ca K line with a spatial resolution of 0." 1 
obtained with the Swedish Solar Telescope (SST) provide evidence that magnetic fibrils, 
originating from network flux concentrations, do span over a large distance above the quiet 
Sun network (see Figs. Ilk. Magnetic field measurements using the He i 10830 A line 



also indicate the presence of a uniform, horizontal magnetic field topology over the inter- 
network at mesogranular scales (Lagg & Merenda 2008). These measurements, presented in 
Fig.[T2] were obtained with the Tenerife Infrared Polarimeter II (TIP-2) mounted behind the 
Vacuum Tower Telescope ( VTT) on Tenerife ( Coll ados et al.|2007] l at a heliocentric angle of 
49° (ji = cos = 0.65). Inversions involving Hanle and Zeeman effect prove the presence 
of a horizontal "canopy" magnetic field on mesogranular scales with strengths of the order 
of 50 to 100 G, similar to the value of the averaged magnetic field of the underlying photo- 
sphere. Both, recent narrow-band Ca K observations of (e.g., |Pietarila et al.|2008 ; Rouppe 
|van der Voort et al.|2 005 ) and magnetic field measurements (e.g., Lagg & Merenda|2008| > 
seem to be in apparent contradiction to Schrijver & Title (2003 ), pointing out the necessity 
of a more detailed analysis on the validity of the concept of the magnetic canopy over quiet 
Sun regions. Nevertheless, the different finding can be fit into a common picture, when tak- 
ing into account the sampled height ranges and a field topology, which is more complex and 
entangled on small scales than usually assumed (see Sect.[5]l. The "classical" canopy might 
be in some ways a too simplified and thus potentially misleading concept. 



4 Numerical simulations of the quiet Sun 

4. 1 Internetwork photosphere 



The solar granulation is now well reproduced by modern radiation (magneto-)hydrodynamical 
simulations. The contrast of continuum intensity or "granulation contrast" is often used for 
comparisons between observations and simulations. For many years, the contrast derived 
from observations were much lower than those found in numerical simulations. One reason 
is the often unknown but crucial effect of an optical instrument and the Earth atmosphere, 
resulting in a significant decrease of the granulation contrast. This problem can partially 
be overcome by using observations with space-borne instruments. Recent observations with 
the Broadband Filter Imager (BFI) of the Solar Optical Telescope (SOT) onboard the Hin- 
ode spacecraft now show higher contrast values. After application of a realistic point spread 
function (Wedemeyer-Bohm 2008, Danilovic et al. 2008 ), state-of-the-art numerical simu- 
lations indeed reproduce important characteristics of "regular" granulation. 

The convective flows in and just above granule interiors advect magnetic field laterally 
towards the intergranular lanes, where the field is concentrated in knots and sheets with up 
to kilo-Gauss field strengths. In the granule interiors, usually only weak field remains, al- 
though in some situations flux concentrations of up to a few hundred Gauss can occur within 



the granules (Steiner et al. 2008). The latter finding is in agreement with the observations by 
|Centeno et al. 1 2007 1 and |Ishikawa et aL (|2008 I. This process of "flux expulsion" has been 
known since early simulations (Galloway & Weiss 1981||Nordlund|1986|). It is now an inte- 
gral part of magnetoconvection simulations (see, e.g .,|Weiss et al.|1996||Stein & Nord lund 



[T998l|Steiner et al.|1998[|Sdia ffenberge r et al.|2005]|y6gler et al.|2005| l. The close-up from 
a simulation by Schaffenberger et al. ( |2005 I in Fig.|l3|illustrates the process. The magnetic 
field in the low photosphere is not only advected laterally. It is also lifted upwards and is con- 
centrated above the reversed granulation layer at a height, which roughly corresponds to the 
classical temperature minimum in semi-empirical models ( [Fontenla et al.|1993| ). There, the 
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Fig. 13 Flux expulsion in a close-up from a MHD simulations by Schaffenberger et al 



(2005 



: Logarithmic 



magnetic field strength in a vertical cross-section (top) and in three horizontal cross-sections (bottom) at 
heights of km, 250 km, and 500 km. The emergent intensity is displayed in the rightmost panel. The arrows 
represent the velocity field in the shown projection planes. The white line in the upper panel marks the height 
of optical depth unity. 



convective overshooting effectively dies out and most of the upward directed flows above 
the granule interiors turn into lateral flows (Fig. [13] In the models by (Wedem eyer et al.] 
2004; Wedemeyer 2003 ), the rms velocity amplitudes are smallest at these heights. The re- 
sult is that the magnetic field is "parked" there and forms a mostly horizontally aligned 
field. It connects to the photospheric flux funnels, which spread out from the intergranular 
lanes below. The enclosed regions below, on the other hand, are virtually field-free with 
field strengths of possibly down to a few Gauss only ( Steiner 2003 ). The field configuration 
around these granular voids was referred to as a dynamic "small-scale canopy" by |Scha£| 
|fenbergeF et al. (2005 2006 ). Virtually field- free granule interiors are very common in their 
simulations. In the more recent simulations by Ste iner et al.|f2~0()iS). this phenomenon is also 
existent but less pronounced, although the field in the granular interiors is still much weaker 
than in the surrounding lanes. The main differences between these simulations is the aver- 
age field strength (10 G and 20 G, resp.) and the injection of horizontal field at the lower 
boundary in the latter simulation. Obviously, the exact occurrence of small-scale canopies 
still depends on details of the simulations and thus needs to be checked by comparison with 
observations. The recent detection of so-called "horizontal inter-network fields" (HIFs) can 
be regarded as observational support for the small-scale field structure seen in the simula- 
tions. It is observed that the horizontal field component in the granular interiors is stronger 
than the vertical component dLites et al.l 19961 lOrozco Suarez et al.|2007l|Lites et al. 12007 
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2008b i. HIFs are also clearly present in simulations ( Schaffenberger et al. 2006; Schii ssler] 



& Vogler 2008 ; Steiner et al. 2008) and are in good agreement with the observations. 



The direction of the horizontal magnetic field, which is continuously lifted to the up- 
per photosphere and lower chromosphere, varies. Consequently, current sheets form where 
different field directions come close to each other. In the simulations by Schaffenberger 
|et al.| (2006) , a complex stacked meshwork of current sheets is generated at heights from 
~ 400 km to ~ 900 km. The lower limit of this range, which is the typical height of the 
small-scale canopies can be considered as the upper boundary of the photosphere. 



4.2 Internetwork chromosphere 

In recent years models have been extended in height to include the chromosphere. Model- 
ing this layer is an intricate problem as many simplifying assumptions, which work fine for 
the lower layers, are not valid for the thinner chromosphere. Rather, time-dependent three- 
dimensional non-equilibrium modeling is mandatory. This is in particular true for the radia- 
tive transfer, for which deviations from the (local) thermodynamic equilibrium should be 
taken into account. Numerically, this is a demanding task. It is unavoidable to make simpli- 
fications and compromises when implementing at least the most important non-equilibrium 
effects in a time-dependent multi-dimensional simulation code. A practicable way is to start 
with simplified models and increase the amount and the accuracy of physical ingredients step 
by step. In their pioneering work, Carlsson & Stein ( 1994, 1995) implemented a detailed ra- 
diative transfer, which was affordable by restricting the simulation to one spatial dimension. 
This simplification made it necessary to implement an artificial piston below the photosphere 
to excite waves as the convection cannot be realistically simulated in one spatial dimen- 
sion. The high computational costs for such detailed radiative transfer calculations forced 
Skar tlien et al.| ( [2000 \ to use a simplified description for their three-dimensional model. 
Nevertheless, their treatment included scattering. Simplifications of the radiative transfer 
are necessary for three-dimensional simulations in order to make them computationally fea- 
sible. This class of 2D/3D numerical simulations cover a small part of the near-surface 
layers and extent vertically from the upper convection zone to the middle chromosphere. 
This way the shock-waves are excited by the simulated convection without any need for an 
artificial driver. The chromospheric layer of these models is usually characterized by intense 
shock wave action, putting high demands on the stability of numerical codes. Wede meyer] 
et al. (2004} made experiments w ith simplified 3D models without magnetic fields, using 



CO 5 BOLD i Freytag et al.|[2002| >. As in the aforementioned simulations, they found that 



overshooting convection in the photosphere triggers acoustic waves that propagate upwards 
and steepen into shock fronts. The result is a dynamic layer above a height of ~ 700km, 
which is composed of hot shock fronts and cool post-shock regions. The gas temperature in 
horizontal cross-sections through the model exhibits highly dynamic mesh-like pattern with 
spatial scales comparable to the granulation. The same can be seen in the follow-up simula- 
tions by Schaffenberger et al. (2005), which include weak magnetic fields (see Fig.|14). The 
gas temperature in the C0 5 BOLD model chromospheres range from about 7000 K down 
to 2000 K, owing to the adiabatic expansion of the post-shock regions. A similar pattern is 
also present in the simulations by Martmez-Syko raet al.| ( p0 08). The temperature range is 
very similar in both models, but the temperature amplitudes differ. Some snapshots of the 
simulation by Martfnez-Syko ra et al.| also show a double-peaked temperature distribution at 
chromospheric heights, but the cool background component is usually much weaker than 
in the C0 5 BOLD model. Possible reasons for the differences are related to the numerical 
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Fig. 14 Horizontal cross-sections through the model by i 
zontal magnetic field component (top), the vertical component (middle row), and the gas temperature (bottom) 
at different heights: z = 0km (granulation), 250km (reversed granulation), 500km, 750km (fluctosphere), 
and 1000 km (from left to right). 



treatment of the radiative transfer in the upper layers. A shock-induced pattern can already 
be perceived in the temperature maps by Ska rtlien et al.| ( |2000) ), although it less pronounced 
due to the relatively coarse grid spacing in this earlier simulation. 

Not only the modeling but also the observation of the shock-dominated layer (hereafter 
referred to as "fluctosphere", see Sect.[5]l is non-trivial. A clear detection in Ca II H, K or the 
infrared lines requires a high spatial, temporal, and spectral resolution, all at the same time. 
A too broad filter wavelength range leads to significant contributions from the photosphere 
below. The fluctospheric pattern is then easily masked by a reversed granulation signal. The 
situation is complicated by the fact that both patterns have very similar spatial scales, i.e. 
roughly granulation scales. This is due to the fact that the generation of both patterns is due 
to processes in the low photosphere. This is illustrated in Fig. [15] which shows preliminary 
synthetic intensity maps in the Call infrared line at A = 854 nm. The maps were calculated 
with the non-LTE radiative transfer code MULTI ( Carlsson 1986) column by column from 
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Fig. 15 Small-scale structure of the solar atmosphere seen in the Ca II infrared line at 854 nm. continuum (left 
column), line wing (middle), and line core (right column). Top row: Synthetic maps based on a simulation 
with non-equilibrium hydrogen ionization; middle row: after application of a PSF and filter transmission; 
bottom: observations with IBIS at the DST (Courtesy of F. Woger). See text for details. 



the model by |Leenaarts & Wedemeyer-B ohm (20061. We use the non-equilibrium electron 
densities, which are output from the time-dependent simulation. The top rightmost panel 
of Fig. [15] shows the mesh-like pattern in the line core, whereas the reversed granulation 
is visible in the line wing (middle column). Even further out in the wing, the granulation 
pattern appears (left column). The mesh-like fluctosphere pattern can be seen Ca H, K, and 
the IR triplet, too. 

A comparison of the line core map with the temperature maps in Fig. [14] shows that 
primarily the hottest regions of the pattern are seen in the Ca intensity, whereas a lot of 
atmospheric fine-structure remains invisible. The hot regions are caused by "collision" of 
neighboring shock fronts, ultimately compressing the gas in the region in-between and ris- 
ing its temperature. This effect enhances in particular the Ca brightness at the vertices of 
the mesh. These small bright areas most likely are observed as Ca grains, while the emis- 
sion along the mesh is so faint that it is hard to detect. The formation of Call grains by 
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propagating shock waves was already explained by (Carlsson & Stein 19971 over a decade 
ago. The fact that their detailed ID simulations closely match observations of grains, clearly 
shows that Call grains are indeed a phenomenon related to shock waves. In ID but also 
in 3D, the formation takes place at heights of ~ 1 Mm above optical depth unity. In both 
cases, the shocks propagate upwards into down-flowing material. The difference, however, 
is that in ID shocks are plane-parallel so that interaction between individual waves is es- 
sentially reduced to shock-merging and shock-overtaking. In 3D, shock wave interaction is 
more complex. And still, the compression zones between shocks - the most likely candidate 
for grain formation in 3D - moves upwards with the waves and thus certainly show very 
similar observational signatures. While it seems to be well established that Ca grains are 
produced by shock waves, some details of the formation process has to be revisited in a 3D 
context. 

However, the grains might just be the "tip of the iceberg". Progress in observational 
techniques and instrumentation now finally allow us to detect the dark details of the fluc- 
tosphere. The middle row of Fig. [13] illustrates this observational effect. A point spread 
function (PSF) has been applied to the synthetic maps. The PSF accounts for a circular, un- 
obstructed aperture of 70 cm diameter and a non-ideal Voigt-like contribution due to instru- 
mental stray-light and atmospheric seeing. Finally, the degraded maps are integrated over 
wavelength with a synthetic transmission filter with a FWHM of 5 pm. The assumptions 
are rather optimistic and represent excellent observational conditions. And yet the resulting 
image degradation has a significant effect on the visible patterns. Obviously, a lower spatial 
or spectral resolution would further suppress the faint mesh-like pattern in the line core. 
Please note that the calculations are still preliminary. A full 3D treatment of the radiative 
transfer and the included scattering, which will soon be possible, might increase the area of 
enhanced brightness. Also it is not clear yet how the possible interaction of the shock waves 
with the overlying "canopy" field would alter the properties of the pattern and its observa- 
tional mesh/grain signature. The resulting pattern nevertheless in many aspects resembles 
the recent observations by F. Woger et al. with (i) the InterferometricBIdimensional Spec- 
trometer ( |Cauzzi et al.|2007|[2008] IBIS) at the Dunn Solar Telescope (DST) of the National 
Solar Observatory at Sacramento Peak (jWedemeyer-Bohm & Wo ger|2008|) and (ii) with th e 
German Vacuum Tower Telescope (VTT) at the Observatorio del Teide I Woger et al.|2 006 1. 
See the lower row of Fig.[l5]for examples of IBIS data. 

Based on the models by Wedemeyer et al. (20 04], we ak magnetic fields were taken into 
account in the simulations by Schaffenberger et al. (2005 , 2006]) and |Steiner et al.| |2008 ) (see 
Fig. |14[l. Different initial magnetic field configurations and strengths from Bq = 10 G to 20 G 
were tried, all resembling quiet Sun internetwork conditions (see Wedemey er-Bohm et al.| 
2005b] for an experiment with Bo = 100 G). The computational domains again comprise 
several granules and extend into the chromosphere, typically to heights of ~ 1400 km. The 
MHD models are very similar to their hydrodynamic precursors with respect to structure and 
dynamics. The ubiquitous shock waves produce a very similar pattern in the gas temperature 
but also shape the small-scale structure of the magnetic field in the upper model atmosphere. 
Consequently, the magnetic field in the fluctosphere is highly dynamic and has a complex 
topology. A look at horizontal cross-sections at different heights in Fig. [14] implies that 
the field in the upper layers is much weaker (\B\ < 50 G) and more homogenous than in 
the photosphere below . On the other hand, the fluctospheric field evolves much faster. The 
horizontal field component Bhor in the range 500 km to 750 km is (i) stronger than the vertical 
one, B z and (ii) has a rather large filling factor there. 

In the small-scale internetwork simulations carried out with C0 5 BOLD, the strongly 
varying surface of plasma B = 1 is found on average at heights of the order of 1000 km 
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to 1400km or even higher, depending on model details. Heights of the same order are also 
found by, e . g . , |Hansteen| p0"07} . The exact location certainly depends on the field strengths 
in the internetwork, which are still under debate. Instead of plasma yS = 1, one can also talk 
about an equivalent surface, where sound speed and Alfven speed are equal. It makes clear 
that these regions are important for the propagation and eventual dissipation. Simulations 
show that this surface indeed separates two domains that differ in their dynamical behavior: 
A slow evolving lower part and a highly dynamic upper part. This is certainly related to 
the finding that wave mode conversion and refraction occurs under the condition of plasma 
P * 1 ( [Rosenthal et"aTp002l |Bogdan et al^[2003) [Caily|p007l |Steiner et al.||2007> . The 
current sheets, which are present below and above the plasma p = 1 surface, differ in their 
orientation. While they are mostly stacked with horizontal orientation in the lower part down 
to the top of the small-scale canopies at the boundary to the photosphere, the thin current 
sheets above plasma p = 1 are formed along shock fronts and can thus show oblique or even 
vertical orientation. 



4.3 Large-scale simulations 

The models described in Sects. |4~Tj and |4.2| do not take into account the large-scale canopy 
fields but rather concentrate on the small spatial scales of quiet Sun internetwork regions. In 
contrast, the simulations discussed in this section comprise larger computational domains. 
To make this possible, one usually has to make compromises such as, e.g., reduce the spa- 
tial resolution or develop efficient numerical methods. Stein et al. (2006) made impressive 
progress by extending the computational box towards supergranulation scales. Their models 
do not include the upper atmosphere but extend deep into the convection zone. Gudi ksen &| 
Nordlund ( 2002 2005a), on the other hand, succeeded in creating time-dependent numerical 
models, which extend from the photosphere all the way into the corona. An important aspect, 
which can be investigated with this kind of models, it the (magnetic) connection between 
the atmospheric layers all the way from the top of the convection to the corona, (see also 
Abbett 2007). Furthermore, extended simulations allow for investigating phenomena that 
are connected to spatial scales between granulation and supergranulation. For instance, the 
simulations by Hansteen & Gudiksen ( 2005 ) and Hansteen et al. ( 2006 ) revealed the forma- 
tion of dynamic chromospheric features similar to dynamic fibrils. Being driven by upward 
propagating waves in the chromosphere, they are an example of the coupling between dif- 
ferent atmospheric layers. Another type of coupling is provided in the form of horizontal 
magnetic flux structures with extensions of a few Mm, which emerge from the upper con- 
vection and rise upwards through the atmosphere. See Martmez-Sykora et al. (2008 ) and 
Che ung et al.|l|2007|> fo r recent examples of flux emergence simulations. The simulations by 



Leenaarts et al.| l 2007[ ) confirm once more (cf . |Carlsson & Stein|2002| l that the ionization 



degree of hydrogen has to be treated in non-equilibrium in the upper atmosphere. Although 
the simulation is two-dimensional, it features weak-field sub-canopy domains with upward 
propagating shock waves and a magnetic-field dominated "canopy" domain above (see their 
Fig. 1). A strong coupling of the individual layers is very obvious. 
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5 An updated picture of the quiet Sun atmosphere 

The results of the previous sections are summarized in a schematic sketch of the quiet Sun 
atmosphere (see Fig. [T6{ | with particular emphasis on the low atmosphere in internetwork 
regions. It is based on (and should be interpreted in comparison with) recent sketches by, 
Judge (2006), and Rutten (2006, 2007) but contains many modifications to incorporate new 
results derived from observations and numerical simulations. 



5.1 The large-scale magnetic field 

The large-scale building blocks of the quiet Sun atmosphere are the magnetic network 
patches, which outline supergranulation cells. The large-scale convective flows (see long ar- 
rows) advect magnetic field to the lanes of the supergranulation. Consequently, the magnetic 
field is highly structured and concentrated close to the "surface" (T500 =1) with kG field 
strengths. The visible result is the so-called magnetic network (see Fig. [7] More recent obser- 
vations with high spatial resolution (e.g., Orozco Suarez et al. 2007 ) reveal that the magnetic 
network patches consist of a conglomerate of smaller magnetic elements or "flux bundles" of 
different field strength with a wealth of substructure. This finding is incorporated in Fig.|16| 
in contrast to earlier sketches that feature the magnetic network as rather massive flux tubes. 
The heights where sound speed and Alfven are equal (c s = c A ), or equivalently where the 
plasma ft = 1, will certainly show large variations, depending on the (local) field strength. 
It may even reach below the surface of optical depth unity (at a reference wavelength of 
500 nm) within strong field concentrations but may stay up at a few hundreds kilometers in 
weaker network patches. And still the variation in field strength and topology, incl. the width 
of the network patches, is even larger than can be presented in the simplified sketch here. 

The magnetic field spreads out in the layers above the patches. Depending on the po- 
larity of neighboring flux concentrations, they can form funnels or connect via loops that 
span the internetwork regions in-between. These two cases are illustrated in Fig. [16] in a 
simplified way. In the classical picture, the large-scale field enclosing the weak-field inter- 
network regions is referred to as "magnetic canopy". The corresponding flux funnels are 
often depicted with a wineglass-like geometry and have their footpoints in the photospheric 
network only. In reality, where the third spatial dimension offers an important additional 
degree of freedom, the field topology is more complex (see, e.g., Gudiksen 2006 ; Peter et al. 
|2006l|Jendersie & Peter|2006[|Schrijver & van Ballegooijen|2005| l. |Schrijver & Title|p003) 
state that as much as half of the field could actually be "rooted" in the internetwork regions. 
From there, it can connect directly to the coronal field or via small loops to the photospheric 
network. The network patches could thus be surrounded by "collars" of loops with spatial 
scales comparable to one or a few granules. Consequently, the concept of a regular canopy 
structure seems questionable. Instead, the field topology should rather be understood as a 
set of individual field lines. Nevertheless, we stick here to the term "canopy" but use it in a 
wider sense. The height of the canopy and the field structure as a whole varies significantly 
from region to region and with time. The height indicators to the left in Fig. [16] should there- 
fore only be used for rough orientation. In principle, the lower boundary of the "canopy" 
field separates two distinct domains: a canopy domain and a subcanopy domain. In reality, 
however, the boundary is certainly less strict than the sketch may imply. Rather, the magnetic 
field of both domains may be interconnected, e.g., by small loops, which extend on granular 
scales (point A). This way, the dynamics of the internetwork photosphere could have a direct 
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influence on the properties of the upper layers, e.g., with respect to wave propagation and 
heating. 



5.2 The canopy domain 

The canopy domain is dominated by (large-scale) magnetic fields. It is this layer, which, 
due to the emission in Ha, appears as a purple-red rim at the beginning and end of a total 
solar eclipse. Therefore, only the canopy domain represents the chromosphere in a strict and 
original sense. At a closer look, a rich fibrilar structure can be seen in chromospheric Ho- 
observations. They are found in rosette-like formations that funnel out from the magnetic 
network below and in many cases connect to neighboring network fields. A few fibrils are 
shown in Fig.[l6]in connection with plasma that is trapped in the chromospheric field. Such 
fibrils and also the larger dynamic fibrils (Hansteen et al. 2006; Langangen et al. 2008a c 
shown at the right in the figure here) are an integral part of the quiet Sun chromosphere 
and even more frequent than can be shown in the 2D sketch here. According to |De Pontieu| 
et al. ( 2007a), fibrils could be the result of chromospheric shock waves that occur when 
convective flows and global oscillations leak into the chromosphere along the field lines of 
magnetic flux concentrations. In general, magnetohydrodynamic waves are an integral and 
ubiquitous part of the canopy domain. (Alfven waves are indicated in Fig.[l6]but represent 
just one of several possible wave modes). Such perturbations can be excited by a number of 
processes, e.g., by the shuffling and braiding of the magnetic footpoints in the photosphere 
by convective flows. As the large-scale magnetic field continues from the lower layers into 
the transition region into the corona above, the whole canopy domain is dynamically cou- 
pled. Again, it must be emphasized that the field topology is certainly more complex than 
can be expressed in the sketch here (see, e.g., Fig.[7]l. Indeed, Sch rijver & van Ballegooijen] 
(2005) state that instead of the plasma-/? surface being closely connected to the (classical) 
canopy, regions with low and high /? can well be mixed up into the corona. 

As already mentioned in Sect. [3] a most obvious constituent of the chromosphere, at 
least when observed at the solar limb, are spicules (see, e.g., De Pontieu et al. 2004). Now, 
two types of spicules are distinguished based on differences in their dynamic behavior ( |De[ 
Ponti eu et al.||2007b"l >. Spicules of type I are the result of shock waves that are excited 
by disturbances in the photosphere (e.g., in connection with p-modes) and propagate from 
there along the magnetic field lines photosphere into the upper layers ( Hans teen et al.|2 006 ; 
Rou ppe van der Voort et al.|2007| l. Spicules of type II, on the other hand, are more dynamic 
but thinner, exhibit higher velocities and have shorter lifetimes (see, e.g., [Langangen et al.| 
2008b ). They are most likely generated by magnetic reconnection events. Alfven waves, 
which by many are considered as an ubiquitous phenomenon in the upper atmosphere, can 
be detected in connection with spicules ( |De Pontieu et al.|200 7d). An example is drawn in 
the upper chromosphere above some vertically orientated current sheets (point C). 

Another ingredient of the sketch are blobs of hot plasma in the corona, although their 
exact position and shape needs further investigation. |De Pontieu e t al. ( 2003 1 showed that 
the emission is not correlated with the centers of flux concentrations. Rather, the emission 
seems to appear at random locations. Although |De Pontieu et al.] refer to "moss" ( |Berger| 
et al. 1999; Fletcher & De Pontieu 1999), which is related to active regions, there is no 
obvious reason why the situation should be different for the quiet Sun corona. Also, hot 
plasma regions like the one marked with "E" in the sketch are certainly not preferentially 
located directly above the middle of an internetwork region. In reality, the entangled and 
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skewed field topology will make such blobs - if existent in the way depicted here - appear 
rather uncorrelated with the field topology of the underlying magnetic network. 

5.3 The sub-canopy domain 

The magnetic field in the sub-canopy domain is mostly weak (see, e.g., |Trujillo Bueno et al.| 
2004j |Orozco Suarez et al.|2007| >, so that the plasma is larger than one in the lower layers. 
There, the field is essentially passively advected by the hydrodynamic flow fields. Con- 
vective motions and overshooting at the "surface" are the fundamental structuring agents, 
making the granulation the dominant spatial scale. Nevertheless, the weak fields in the sub- 
canopy domain most likely connect at least partially with the stronger canopy field. This 
feature is taken into account as integral part of the atmosphere sketch. Unfortunately, the 
presentation remains rather speculative at this point as many details of how and where the 
connections exactly take place are still unknown. 

Beside the magnetic field, the consequences of convective overshooting allow to divide 
the subcanopy domain into layers with distinct dynamics (from bottom to top): low photo- 
sphere, middle photosphere, high photosphere, fluctosphere. 

The lower and middle photosphere exhibit the visible imprints of the solar surface convec- 
tion. The granulation in the low photosphere is directly produced by small-scale convection 
cells (see, e.g., Nordlund & Dravins 19901, while the reversed granulation in the middle 
photosphere is a second-order effect. Gas is brought up by convective overshooting in the 
granule interiors, adiabatically expanding and cooling. It streams down again in the inter- 
granular lanes, where it is compressed and heated. In addition, p-modes, i.e. global oscil- 
lations, and local acoustic events are important ingredients of the photospheric dynamics. 
Recently, Stra us et"aT] ([2008 ) presented new support for the idea that gravity waves could 
play an important role, too. 

The usually weak magnetic field is brought up from the convection zone below and/or 
possibly locally generated by small-scale dynamo action close to the surface. In the pho- 
tosphere, the weak field is more or less passively advected towards the intergranular lanes 
but also towards the upper photosphere. The resulting field concentrations in the lanes be- 
come visible as very small and confined structures, e.g., in G-band images (see, e.g., [de] 
Wijn et al. 2008 for a recent example). In general, the internetwork field in the photosphere 
exhibits significant inclination and mixed polarity (see, e.g., |Martmez Gonzalez et aTp 008 ; 
|Orozco Su arez et al. 2007). The granule interiors may become virtually field-free if there 
is no supply of magnetic fields with the warm convective upflows. Such voids are enclosed 
by small-scale canopies. Over most of the granulation, the horizontal field component is 
stronger than the vertical. This effect is observed as "horizontal internetwork fields" (HIFs). 
Magnetic field can emerge also in the form of small loops, which may have footpoints even 
within a granule (see point B in Fig.|16[l. This process, which was observed by |Centeno et al.| 
( 2007J, most likely adds to the accumulation of field above granules. In addition to emerging 
loops, [STein & Nordlund|p006| l report on flux that is submerging and thus disappears from 
the surface. 

The upper photosphere marks the boundary between the photosphere, which is controlled 
by the effects related to convective overshooting, and the wave-dominated layer above. This 
boundary can roughly be placed at the height of the classical temperature minimum. There, 
the temperature structure appears smoothed out and less structured than above and below; it 
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is here that the average temperature amplitudes are smallest. It is roughly the height where 
the UV continuum at 160 nm is formed (cf. TRACE passbands). The upper photosphere 
is the layer, where the small-scale canopies have their top and where stacked (horizontal) 
current sheet become most obvious. This layer can be seen as a kind of (dynamical) insula- 
tion between the internetwork photosphere and fluctosphere. This effect becomes obvious in 
simulations when starting from an initial condition which feeds in field at the lower domain 
boundary. The photospheric field is built-up rather quickly but the field above only after a 
time delay because it only slowly spreads into the strongly subadiabatic stratification of the 
upper photosphere. 



The fluctosphere: The shock-dominated domain in subcanopy internetwork regions (see 
Sect. |4.2j>, is referred to as "fluctosphere" by |Wedemeyer-B6hm & Woger| l |2008| >, while 
Rutten (20071 uses the term "clapotisphere". It is located between the photosphere and the 
part of the chromosphere visible in Ho-. It is composed of propagating and interacting shock 
waves (with weak field only) and intermediate cool post-shock regions. Ideally, the wave 
fronts would expand spherically, while moving in vertical direction. In reality, they are de- 
formed by running into an inhomogeneous medium of downflowing gas, which was shaped 
by precursory wave trains. The horizontal expansion of the fronts inevitably causes interac- 
tion between them. A visible result is the formation of Ca grains at heights, which tradition- 
ally would be assigned to the low chromosphere. The waves are excited in the photosphere 
below via different processes, which are related to convection (e.g., exploding granules), 
overshooting, and p-modes. The magnetic field in the fluctosphere is rather weak and is 
therefore mostly passively shuffled around by the shock waves. The result is a very dynamic 
and entangled field. The strongly varying surface of plasma yS = 1 or in this context better 
c s = ca is most likely located at heights of the order 1000 km to 1500 km or even higher (see 
Sect. |4.2| >. There, the conditions allow wave mode conversion, so the parts of the fluctosphere 
below and above can show a somewhat different dynamical behavior. In the upper part, the 
weak fields become more important and rapidly moving filaments of enhanced field strength 
are generated. The propagating shock waves nevertheless remain the dominating structuring 
agent. A consequence, however, is that the current sheets are only stacked at plasma f3 > 1 . 
Above, they are less regular as they are formed in the narrow collision zones of shocks, 
where the magnetic field is occasionally compressed. This shock-induced magnetic field 
compression might qualify as a (minor) heating process with potential consequences for the 
chromospheric energy balance. 

The fluctosphere is not directly visible in Ho- (in the line core at least) and is thus not 
a part of the chromosphere in a strict sense. It seems advisable to reserve the term chromo- 
sphere for the fibrilar canopy domain as visible in Ho- (or in the very line cores of the Ca II 
lines). However, the fluctospheric shock waves could still leave an imprint in chromospheric 
diagnostics by interacting and penetrating the canopy field. On the other hand, the flucto- 
sphere is also no part of the photosphere, although causally connected via the shock waves 
that propagate upwards from the low photosphere. The fluctosphere could be regarded as a 
second-order effect only, in contrast to the granulation and reversed granulation, which are 
direct consequences of the solar surface convection. 



5.4 Shock waves meet the "canopy" 

Some details of Fig.[l6]concern the interplay of propagating waves and the magnetic canopy. 
There are certain zones in these magnetic structures that act as mode conversion zone (Bog- 
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|dan et aL|2003[ [Ca lly 2007), e.g., converting incoming acoustic waves into other modes, 
such as fast and slow magnetoacoustic waves. It is thus possible that such converted waves 
continue to propagate along the canopy field lines as some kind of "canopy waves". For sim- 
plicity, such a zone is marked by "wave conversion" at the "outer" boundary between canopy 
and subcanopy domain in the figure. Generally, such zones can be located everywhere in the 
structure where sound speed and Alfven speed are of equal magnitude. 

Furthermore, refraction and even reflection of waves can occur in such zones. As for 
the mode conversion, details depend on the relative orientation of the incoming wave and 
the magnetic field (Hasan et al. 2008]l. A wave can remain barely affected by the field when 
traveling perpendicular to the field lines, e.g., upwards in a vertical flux concentration. On 
the other hand, significant (relative) inclination can even result in total internal reflection for 
some wave modes (Rosenthal et al. 2002). In general, it can be assumed that the (acoustic) 
shock waves coming from the fluctosphere are guided by the magnetic canopy (e.g., point D 
in Fig.|16|l. Consequently, waves might follow the canopy field upwards and compress and 
heat the gas trapped between chromospheric "funnels" (point E). In closed loop regions, 
strong waves could push into the canopy from below and compress the magnetic field (loca- 
tion F). Depending on the local field configuration and the properties of the incoming wave, 
such an event could eventually trigger reconnection events. It could contribute to chromo- 
spheric heating. It certainly would not be only limited to the locations indicated in the sketch 
but occur more often in complex 3D field configurations. On the other hand, a regular closed 
structure as in the figure could possibly refract the waves from below such that they are "fo- 
cussed" in the top of the subcanopy domain, amplifying their effect on the canopy field. A 
possible - although speculative - result could be the triggering of "nanoflares", although 
they are initiated by other mechanisms at other locations, too. That the upwards propagating 
waves interact with the canopy field is implied by observations in Ho-. The dynamic behav- 
ior of the chromospheric fibrils is reminiscent of strings that sway back and forth in reaction 
to the quasi-continuous impact of waves from below (point F). Under certain conditions, the 
shock waves might actively deform the field configuration of the magnetic canopy. Magne- 
toacoustic waves can already enter network flux concentrations in the photosphere, where 
the inclined magnetic field lines act as "magnetoacoustic portals" (Jefferies et al. 2006). 
The observations of so-called "acoustic shadows" provide observational evidence for the 
interaction of acoustic waves with the field around network footpoints (Krijger et al.|2001"[ 
|McIntosh et al.|2001[|McIntosh & Judge|2001| >. 

5.5 Probing the upper atmosphere 

With the currently available diagnostics for the chromosphere, observations of the sub- 
canopy domains are problematic. The Ha line core samples only the "canopy domain", 
whereas observations in the line wing reveal a background that most likely is dominated by 
the reversed granulation at much lower heights. It seems questionable if the layer in-between 
- the fluctosphere - can be observed in the Ha line wing at all in internetwork regions. Po- 
larimetric measurements in the He 10830 A line (see Sect. [3} principally allow for the de- 
termination of the magnetic structure in a slab located between 1000 km and 2000 km. The 
formation of this line requires coronal illumination in the UV, resulting in complete absence 
of any photospheric contamination. However, the main contribution in the He 10830 A line 
comes from layers slightly above the fluctosphere. The Ca H & K and IR lines in principle 
would allow observations of the fluctosphere if very narrow filters are used. Otherwise, the 
detected intensity is "contaminated" with radiation from layers below. Very often, Ca ob- 
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servations with too broad filter prominently show the reversed granulation (see Figs.jTJl and 
[TTj i, which is easily mistaken as chromospheric signal. Very narrow filters, on the other hand, 
make it necessary to properly correct for Doppler shifts. A solution is fast scans through the 
Ca II IR lines with new imaging polarimeters such as IBIS ( Cauzzi et al. 2007 , 2008 , Kleint 
et al.|2008|l or CRIS P ( |Scharmer et al.|2008fr , or spectro-polarimeters like SPINOR l |Socas 



Nav arro et al.|2 006 ). The extended formation height ranges and the non-equilibrium condi- 



tions, under which the inner parts of these lines are formed, complicate the interpretation and 
the derivation of the atmospheric structure. A promising alternative are the (sub-)millimeter 
continua, which will become accessible with the Atacama Large Millimeter Array (ALMA) 
a few years from now. Although technical details of this new type of observation render the 
construction of brightness temperature maps a certainly very complicated task, the scien- 
tific results could significantly contribute to our understanding of the solar atmosphere at 
chromospheric heights (Loukit cheva et al.] 2008 , Wedemeyer- Bohm et al.|2007| . 



6 Conclusions 

The solar atmosphere is a very dynamic and inhomogeneous multi-scale system. Its individ- 
ual components are coupled; some of them even show a kind of hierarchical self-similarity. 
Examples are the observational imprints of sub-surface convection, with a continuous spec- 
trum of scales from below granulation scales to above supergranulation scales, and magnetic 
fields, which also exhibit similar features over a large range of spatial scales. 

Despite great progress on the theoretical and observational sides, which go hand in hand, 
we are still missing an ultimate, comprehensive picture of the quiet Sun atmosphere. But at 
least we can now see what is needed for a corresponding numerical simulation. First, the 
computational domain should be large enough to encompass a few supergranulation cells 
while the spatial resolution must still be high enough to capture important processes that oc- 
cur on scales smaller than granulation. The vertical couplings make it necessary to consider 
an extensive height range. The corona and chromosphere can only be treated realistically 
when including the important driving motions in the layers below, i.e. in the photosphere and 
(at least) the upper part of the convection zone. While many simplifying assumptions can be 
made for the lower parts of such a model, the layers above the (middle) photosphere require 
a numerically complicated and thus computationally expensive non-equilibrium modeling 
approach, e.g., a realistic treatment of hydrogen ionization etc.. The production of such a 
comprehensive model - and analogous models for, e.g., active regions - is thus very involved 
and can be regarded as one of the current challenges in (computational) solar physics. 

On the observational side, we must continue to push forward the instrumental possi- 
bilities towards higher resolution in the spatial, temporal, and spectral domains - all at the 
same time. In addition, the (further) development and exploitation of advanced diagnostics 
is needed to derive a seamless tomography of the atmosphere as an integral phenomenon. 

Until we succeed to reach these ambitious goals, we are left with a number of open 
questions. Of particular interest for the quiet Sun are, amongst others: 

- How does the weak internetwork field connect with the stronger network field? What 
does the magnetic field look like just below the "canopy"? And can we talk about a 
"canopy" even in a wider sense after all? 

- How do the propagating fluctospheric shock waves interact with the stronger ("canopy") 
field? Is it mostly a "passive" refraction / reflection at the "boundaries" of flux con- 
centrations (e.g. [Rosenthal et al.| |2002; Stei ner et al.||2007| or also "active" distor- 
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tion/displacement/compression of magnetic field? How and where does mode conver- 
sion take exactly place under realistic conditions? 
- The question of the coupling between the atmospheric layers is closely connected to the 
heating mechanism question or, better said, to the question of the atmospheric energy 
balance, not only for the Sun but also for other stars. Amongst other things, the ongoing 
controversy concerning the heating mechanism of the quiet Sun chromosphere (e.g., 
|Fossum & Carls son 2005 ) has important implications for stellar activity in general. 
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